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The Quantum Stochastic Differential Equation Is Unitarily Equivalent 
to a Symmetric Boundary Value Problem for the Schrodinger Equation 
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Abstract. We prove that the solution of the Hudson-Parthasarathy quantum stochastic differential equation in 
the Fock space coincides with the solution of a symmetric boundary value problem for the Schrodinger equation 
in the interaction representation generated by the energy operator of the environment. The boundary conditions 
describe the jumps in the phase and the amplitude of the Fourier transforms of the Fock vector components as 
any of its arguments changes the sign. The corresponding Markov evolution equation (the Lindblad equation 
or the "master equation") is derived from the boundary value problem for the Schrodinger equation. 
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§1. Introduction 

j ' Our interest in the theory of quantum stochastic differential equations is primarily due to the fact that 
it can be used to represent solutions of the Markov evolution equation, which generalizes the Heisenberg, 
\£) Kolmogorov-Feller, and heat equations [1-4]. 

O Solutions of Markov evolution equations are studied in the theory of open quantum systems [5-8]; 
they are called quantum dynamical semigroups. One of the consequences of the stochastic representation 
of quantum dynamical semigroups in von Neumann algebras is a simple relationship between the con- 

j» servativity of a quantum dynamical semigroup and the isometric property of a solution of the quantum 

gtochastic differential equation. The conservativity (or the unital property) of a quantum dynamical semi- 
""^3 group means that the unit of the operator algebra is preserved; this corresponds to the preservation of the 
total probability by the solutions of the equations for the transition probabilities of the stochastic processes 
j 1 , in the classical and quantum cases. Violations of this property are related to critical phenomena such as 
O explosions, escape of the solution to infinity, accumulation of infinitely many discontinuities, nonunique 
S solvability of the Cauchy problem, the existence of growing solutions to formally dissipative equations, 
r ^ etc. [9-10]. The study of conditions necessary and sufficient for the existence or nonexistence of such 
L* phenomena is currently of considerable interest [11-18]. 
. — Sufficient criteria for solutions of quantum stochastic differential equations to be isometric or unitary 
were obtained in recent years mainly by perturbation methods [19-21]. However, symmetric operators 
^ Responsible for these properties of solutions were not found. The difficulties encountered in this approach 
include the violation of the group property by solutions of quantum stochastic differential equations [22] 
and the fact that their formal generators fail to be symmetric on domains consisting of smooth functions. 
More precisely [23], the formal generators for the Schrodinger equations unitarily equivalent to quantum 
stochastic differential equations have the form of dissipative operators perturbed by singular quadratic 
forms in the sense of [24-25] . Since the dissipative part of the generator has an unbounded anti-Hermitian 
part, the general theory of self-adjoint extensions of symmetric operators bounded below and perturbed 
by singular quadratic forms [24-26] does not apply directly to this case. 

In the present paper, we consider quantum stochastic differential equations whose solutions in the in- 
teraction representation are limits of solutions of a family of Schrodinger equations in the Fock space 
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[23]. The fundamental property of functions belonging to the range of the resolvent of the limit uni- 
tary group is that they experience amplitude and phase jumps at the points where the quadratic form of 
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the generator has singularities. Evaluating the action of the generator of the limit group on functions 
with amplitude and phase jumps, we obtain a symmetric operator; its self-adjointness can be verified 
independently under assumptions less restrictive than those used to construct the resolvent explicitly. 
Thus, we prove the equivalence of the quantum stochastic differential equation and the boundary value 
problem for the Schrodinger equation in the Fock space. Straightforward computations permit one to derive 
the Markov evolution equation directly from the boundary value problem for the Schrodinger equation. 

We start from a simple example showing the distinction between the formal weak limit, the strong 
resolvent limit of a sequence of generators of unitary groups, and the action of the generator of the limit 
unitary group on the range of its resolvent. 

§2. The weak limit and the resolvent limit of a 
family of generators of one-parameter unitary groups 

Consider the one-parameter family of unitary groups exp{itH a } = in L 2 (M.) given by 

U^ a) Tp{x) = V>(x-f)expjzA J drV a {x -t + r)^, i,Ael, 

where V a (x) = exp{— x 2 /2a} , a e M + . Obviously, V a (x) — > 5{x) as a — > +0. Therefore, 

the weak limit of the family of essentially self-adjoint operators H a = id x + XV a (x) is described by the 
bilinear form 

H w [<p, V] - (<p, H w il>) = i(<p, V') + A^(0)V(0), (2.1) 
which is well defined on W • O n t ne other hand, 

-t 

drV a (x -t + r) -> 7 [0 . t) (x), 

o 



where It{x) is the characteristic function of a Borel set T C K. Hence, the strong limit of the family 

U[ a) is 

lim U [ t a) ^{x) = UMx) = e ltH «V(x) = Mx - t)e lX ^^ [x) . 

Note the identity e*^ '*)^ = (e lX — l)I[o,t)( x ) + 1- Therefore, the bilinear form of the resolvent limit 
Jin = r-\\mH a is well defined on W 2 1 (M) and is given by 

H R [<p,1>] =r 1 lim^( V ,^V) _ = i(<p,i/>')+i(e iX -l)<p(0)il>(0). (2-2) 

i *U CxXj t — 

Comparing (2.1) with (2.2), we conclude that H_r = r-\\m.H a ^ w-Ym\H a = H w . 

The range of the resolvent of the limit unitary group Ut is the natural domain of the generator Hr 
and can be described explicitly in our example: 

POO POO 

R„ip(x)= dte-^7p(x-t) + e(x)(e iX -^e'^ dte~^(-t), 



where 8(x) is the characteristic function of the semiaxis M + . This structure of the resolvent means that 
the functions belonging to the domain of H# have a phase jump at the origin x = 0: 

lim Rutp(x) = e lX lim R^{x). 

x^+0 x^ — Q 

Therefore, the domain T>\ of the generator H« consists of the functions 

V>: tp G W 2 1 (R \ {0}), lim ip(x) = e iX lim <p(x). (2.3) 

x^+0 x^ — 
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The operator acts as id x for x ^ . For functions tp e V x , the left and the right limits at the origin 
exist by virtue of the embedding W2 (K\ {0}) C C(R\ {0}) . The symmetry property of Hr follows ^from 
integration by parts and from the identity 

<p(x)ip(x)\~l = 0, 

which holds for ip , tp e T>\ ; the sclf-adjointness follows [27] from the fact that for p = ±1 the equation 
(Hfl + ip,)ip(x) — f(x) , x 7^ 0, with the boundary condition (2.3) is solvable in T>\ for every right-hand 
side / e L 2 (R). 

§3. The weak and the strong resolvent limit of 
solutions of the Schrodinger equation in the Fock space 

Suppose that Ti is a Hilbert space, r^Z^M)) is the symmetric Fock space and h = H <g> r s (_L2(K0) 
is their tensor product. For v,g <G L 2 {M.) , by A + (v) and A(g) we denote the standard creation and 
annihilation operators in r s (L 2 (K0)- Consider the family of Schrodinger equations d t ipt = iHip t with 
self-adjoint Hamiltonian H = H ®I + I®E + H lnt in h: 

H int = K® A+(g)A(g)+R®A+(f) + R*®A(f), E = j ' uoa+(uu)a{u) dcu. (3.1) 

^From now on, it is assumed for simplicity that the operators Hq, K, and R € C(H) commute and have 
a joint spectral family E\ . More precisely, let 



H = Ju(X)dE x , K = JxdE x , R = Jp(\y*^ 



dE x , 



where v, p, and $ are measurable real functions corresponding to the operators H 0l K , and i?, so that 
the operators H and K are self-adjoint and R is normal. 

We denote by Pt(X) the one-parameter group of unitary operators in L2QR) with the generator iV(A) = 
ui + X\g)(g\ . The unitary group U t = cxp{iHt} can be constructed by transforming the Hamiltonian to 
the canonical form [28] , or by using the interaction representation generated by the operator J ® E + K ® 
A + (g)A{g) [23]. Theorem 3.1 describes the action of U t on coherent vectors. 

Theorem 3.1. The unitary one-parameter group Ut = cxp{itH} , where H = Ho<E)I + I<£>~E + H lnt 
and H lnt is defined by Eq. (3.1), acts as follows: 

U t h ® i,{v) = J e-W* dE x h ® ip(y t (X)) exp{ip(A)e-^ A ) J (/, v s (X)) d*| , (3.2) 
where ^ 

v t (x) = P t {x)v + ip(xy*^ [ P s (x)fds. 



Let L^~(R) C W^R) be the set of functions with positive absolutely integrable Fourier transform; 



W.J W is the Fourier transform of the Sobolev space W r 2 1 (IR) . Let /, g e L+~(R) be real functions such 
that /(0) = 5 (0) = 1/v^F, and let f( a )(u) = f(atu) , g( a \u) = g(atu) . 

Lemma 3.1 (about four limits). Let P^ a \x) be the one-parameter unitary group in L2OR) with gen- 
erator N a (X) = uj + X\g^)(g^\ , and let ttt = Ft->ulT(t)FZ,->t be a family of projections in H(L 2 { 
Then the following limits exist as a — > : 



1) lim/ d8tf a \pM(\)fW) =(2-iX)~ 1 
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2) w-lim /" t dflPi a )(A)/W(w) = e <w *7 [0 , t )(a;)(l-iA/2)- 1 ; 

Jo 

3) w-\im(g( a \P t {a \x)v) = (1 - iX/2)-^*^ t v; 

4) s-limP t (Q) (A) = exp{iZ(A)Tf ((M) } = P t (A) , exp{iZ(X)} = (2 + iA)/(2 - *A) . 

Let us consider the family H a of self-adjoint Hamiltonians (3.1) parametrized by the functional argu- 
ments f( a \u>) , g^ a \u>) of the creation and annihilation operators. Lemma 3.1 in [23] justifies the passage 
to the limit in Eq. (3.2) as a — > +0. The limits l)-4) correspond to the passage to the limit in the four 
components of the solution (3.2) depending on a. By substituting l)-4) into Eq. (3.2), we obtain the 
limit unitary group XJ% = cxpjzHi} = s-lim a _»o exp{iH a i}: 

U t h ® ip{v) = J e iH W* dE x hi>(e iZ ^ a ^e ilJJt v + j^p{X)e^I [Q ^ 

x ex P {p(A)e-^W^ x (7 [ o, 4 ),e i ^)}, (3.3) 

where H(X) = v{X)+ip 2 / (2 — iX) . Let W and L be the operators with spectral densities (2 + zA)/(2 — iX) 
and 2i/(2 — i A) p(X)e 1 ^^ , respectively. Then 2ip(X)e~ l< ^^ / (2 — iX) is the spectral density of the operator 
—L*W, where W = (2 + iK) /(2 — iK) is the Cayley transformation of the self-adjoint operator 2K , and 
L = 2i/(2 — iK) R is a densely defined operator such that doml D domi?, and H(X) is the spectral 
function of the operator 

iG = H - \l*KL + l -L*L. 

In the following, we suppose that the operator — G is the generator of a strongly continuous one-parameter 
contraction semigroup W t = cxp{— Gt} in H, and moreover, 

D = domff ndomi'L C domG C domi, G*tp + Gtp = L* Lip VtpeD, 

where the operator H = —Hq + jL*KL is symmetric on D . In this notation, the bilinear form 

H[g ® V(/) , h ® xp(v)\ - t lim o \j t (9® V>(/) . U t h ® V(«)) , 

where g , h e D C H and f ,v £ (R) i acts as follows: 

H[<? ® V(/),/i® V>(«)] = e (/ '" )t2 Gft) w - i(g,Lh) n J{Q)+i(L9, Wh) n v(0) 

+ {g,h) n j dw f(u>)g(u>)u> + i(g,(I- W)h) n W)v(0?j ■ (3-4) 

The bilinear form (3.4) has the regular dissipative component ^L*L®I and the singular component [24- 
25] 

—iL®A + (— ^=) + iL'W ® a(— ^=) — W) ® A + (— ^=)a(— ^=) , 
vanishing on the subset 2? 

V = {$:$ = ft ® V(«) , V G W^QR) , «(0) = 0,fieD}, 

which is total in h. At the same time, the weak limit of the sequence H a is described by a bilinear form 
with different regular and singular components: 

\imH a \g®^{f),h®^{v)]=e^' v ^(^g,H h)n + {9,h)n j codujJ(w)g(w) 

+ (g, Kh)nW)v(0) + (.9, Rh)nW) ~ (Rg , h) n v(0) 
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Just as in the example in §2, here we have 

r — lim H, v ^ w — lim H Q . 

Let us consider the set function u(s,t) — J s Ut- s Jt i where Jt is the unitary one-parameter group 
generated by the operator E = J dui uoa + (uu)a(uj) [7], which acts in V s as Jtip{v) = ip(e lu ' t v) . Taking into 
account the fact that J^tp(v) = %l)(e~ wt v) , ^from Eq. (3.3) we obtain 

u(T)h ® ip(v) — f e iH ^ mesT dE x h®ip(e iz ^ T v + ip{X)e^ W — 2 — T T ) 

xc^{zp(X)e-^ J ^(I T ,v)}. (3.5) 



For disjoint sets, the family of commuting operators u(T) satisfies the composition law u(T\ U T%) = 
u{T\)u(T2) , and the differential of the bilinear form (g®ip(f), u(0, t)h®tp(v)) satisfies the weak quantum 
stochastic differential equation [2] 

d(h®ip(v),u(0,t)h®ip(v)) =i(h®ip(v),u(Q,t)H(dt+)h®ip(v)), (3.6) 

where 

iH{T) = M{T) = J dtJ t (U-E®I)J*, 
M(T) = -G ® mcs T + L ® A+{T) - L*W ® A{T) + (W — I)® A(T). 

Thus, the following statement holds. 

Theorem 3.2. The family of solutions of the Schrddinger equation with Hamiltonian (3.1) in h strongly 
converges, up to the unitary transformation J t * , to the solution u(0,t) of the stochastic equation (3.6), 
u(0, t) — s -lim a ^Qlj[ a \jl . 

In what follows, we shall see that, quite unexpectedly, the natural domain of the generator of the 
group U t (the range of the resolvent) does not contain functions on which the bilinear form (3.4) is well 
defined: the Fourier transforms of the Fock components of the resolvent experience amplitude and phase 
jumps at the points where their arguments change sign. 

§4. Surprising properties of the resolvent 

Consider the following vector $ e h = H ® V s (L 2 ( 



<S> = R fl h®i>(v)= dte-' a Uth®il)(v) = {* n (u)}, $„(■): R n -> H, u = {u u . . . , u n }, 
Jo 

with components (3.3) 

$ n (w) = J dt exp|-/i(G + t - L*W u(-T)dr|^ n , t (w), 

n 

Vn.tM = ]J((W - l)7r [0;t) e^*w(a; fe ) + e iUkt v(w k ) + LI [0it) (u k ))h, 
i 

where L , W , and G are the above-described commuting operators with spectral densities 

L(X) = 2ip(X)e- i ^ x \2-iX)-\ W(X) = e iZ ™ , G(X) = -iv{X) + . 
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Let (f n ,t(^) be the Fourier transform of <p n ,t{w)i 

n 

<fn,t{ T ) = X\{{ w - I)I[0-,t)( T k)v(Tk - t) +v(r k -t) + LI [0tt) (T k ))h, 
1 

where r = {n , . . . , r„} . Let Kc{l,...,n}, and let /C a be the complement of /C . Set 

P SW = II ^ w - 7 ^ Tfe - *) + L ) 7 [o,*)(^) e B(h). 



feG/C 

Then 



; »,«w=Eteiw n (4.i) 

Obviously, the functions Pj^\ (r) have discontinuities of the first kind at the points where the variables 
Tfc change the sign: 



Jim o PW(r)=/ X ;.(fc)pW(r), 

lirr iPg(r) = /^(fc)Pg(r) + ((W - I)ff(-t) + i)P^j it (r). 



(4.2) 



Therefore, Eq. (4.2) implies 

= -((W- m-t) + L)P^ t {r)I K {k). (4.3) 
Let us calculate the jump of ^ n ,t(r) at the points where r k changes sign. Note that 

lim (p n ,t{ T ) = «(-i)^™-i,t( r i> ■•■ , Tfc-i, Tfc+i, ... ,t„). 

Taking into account Eqs. (4.1) and (4.3), we find the amplitude and phase jumps of the functions from 
the domain of the infinitesimal operator of Ut : 

lim !p nt ( T )=W lim ^„, t (r) + £p n -i,t(n , . . . , r fe _i , r fe+ i , . . . , r n ). (4.4) 

Let XV, l = -D® L S (W 2 1 (M\ {0})) C h be the vector subspace of elements satisfying condition (4.4), and 
let A(5±) , A(5±) , and N be the operators acting on the Fock vectors by the rule 

I " 

(*,A (5±W = £ hm L 0}) ^ n <M*»,*»)< =±e , 

1 fe= i"'(K\{0}) m ^ k 



fe=i 



In this notation, the boundary condition (4.4) in F s acquires the form 

(N + l- 1 (I ® A(<5+) - W ® i4(J_))* = L <g> I . (4.5) 

This condition is imposed on the Fock vector components ^2, • ■ • and extends condition (2.2) to the 
Fock space. 
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Let us prove that the operator 

H = iG®I + I®E + iL*W <g> A(5-), E= [ dr a + (r)a(r)^ r , (4.6) 



/R\{0} 

is symmetric in T>w,l ■ Let $, G T>w,l , and let B be a Hermitian operator such that domB ® 7 D 
T>w,l ■ Integration by parts yields the following identity, where the integrated terms are expressed via the 
operators A(6±): 

($, B <g> m) - (B <g> E$, *) = B <g> (A(5_) - A(J+))#) 

oo n 

— < )/! ^ — / - 



oo 1 n 

1 n - fc=l- y ( R \{°»"- 1 m^fe 



Using the boundary condition (4.5) for and we find the integrated terms in (4.7): 

»($, S ® (A(<y+) - A(<y_))*) = (TU*SVK - B)A(J) - *) 

BL^>) + I(WA(6-)&, BL^>) + i(L<f>, BWA(6-)V>). (4.8) 

In particular, Eq. (4.8) is simplified for B = I as follows: 

/(*, / ® (A(J+) - A(<y_))*) = »(*, L*Ltf ) - (iL*W ® A(<y_)*, *) + (*, ® ). 

Since iG — iL*L = (iG)* , we obtain the following identity, which means that the operator H is 
symmetric in T>w 

= ((I® 25$,*) + ($,{iG<8> J + *L*W®j4(J-)}tf) 

-t($, J® (A(5+) - A(<y_))*) 

= ((/ ® 25$, tf) + ($, (iG)* ® Jtf) + (*L*W ® A((5_)$, tf) = (/?$, *). 

Let us find how the generator of the group U t acts on the range of the resolvent. Set <J/ <G h and 
$ = R^h ® V ; (' y ) ■ By the definition of the generator, we have 

^ 1 J 1 /"OO 00 1 /" 

(*,fl"$)= lim - — (*, t/,$) = - / dt^T — dr 
s^+o ids i Jq ^ n\ 7(R\ {0 })" 

V as / h s=o 

Note that the functions tp n ,t(T~) depend on the differences Tfc — t. Therefore, 

d _ "a 

T^.iW =-^^— ¥>„, t (T) =t£7^„, t (T). (4.10) 

fe=i Tfe 

On the other hand, 

A(5_)^ n , t (T) = nw(-t)^ n _i lt (r) (4.11) 

by the definition of A(5-) . Now, taking into account definition (4.9) and identities (4.10) and (4.11), we 
obtain 

/■oo / \ r oo °° -i /> 

difVciGe-^+^-^^/o^^/j +/ dtY^-l dr 
Jo V Jh h ^ «' J(E\{o})" 



<9r fe 



<Pn,t{ T ) 



x C^„(t) , e-CG+Mjt-ii'w/o' «(-■>-) <Jt ^G + »L*W«(-i) + 7^ ■ 

v v fe=1 

= (f,{iG + ii*IU®^((5_)+/®i?}$) h ; 
that is, the generator 27 of the group U t has the form (4.6). Thus, we have proved the following theorem. 
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Theorem 4.1. The symmetric operator 

H = iG®I + I®E + iL*W® A(6- ) 

in T>w,l, where G = iH + jL*L, H = jL*KL — H , is the generator of the one-parameter unitary 
group U t ■ 

We must point out that the verification of the symmetry property did not rely on the assumption that 
L, G , and W commute and can readily be extended to operators of the form 

H = iG + I®E + i'%2LtW l , m ® A m (6-) (4.12) 

l,m 

with the boundary condition 

(N + l)- 1 (i ® MS+) - w i,m ® A m (6-)\ * = U ® I* , (4.13) 

where T4 7 = {W; ;TO } is an M x M-matrix with entries in B(H) such that W*W = I, and {A m (g) : g e 
Z^QR); 1 < m < M} are the annihilation operators in r s (L2(R M )) , which commute for different L 

§5. The Markov evolution equation 

In conclusion, we describe a derivation of the Markov evolution equation from the Cauchy problem for 
the Schrodinger equation 



d_ 

It 



*(*) = (-G + U ® E - Y, LW,m ® ^ m (<5-)) *(t) 



with the boundary conditions (4.13). Let B e B(H) be a Hermitian operator, and let h, g e I? . Consider 
the equation 

(G, P t (B)h) n = (U t g ® *(0)|B ® /|C/ t /i ® *(0)) h 
for the mean value. From (4.12) we have 



(g, P t (B)h) H = -((g + J2 LfWi^AniS-)) U t g ® *(0)|B ® J|£/*/i ® *(0) 



I'm 



- (u t g ® *(0)|B ® /| (g + ^ i^W^Ar^-)) t/ t /i ® 4-(0)] 

+ ^ tfl ®*(0)|B® (A(J+)-A(J_))|t^/i®*(0))^. (5.1) 



Now we can use identity (4.8), which in this case can be rewritten as follows: 



+ ]T (£;<&, Bi ; *) + ^ ((Wi.m^OJ-)*. BL^) + (L,$, Bl^, mJ 4 m (5_)*) N ) . (5.2) 
i i,m^ ' 



518 



and A(£_)*(0) = 0, ;,£rom Eqs. (5.1) and (5.2) wc 

-(Gg, Bh) - (g, BGh) + J2( L i9, BL l9 ). 

i 

Thus, we have obtained the generator £(■) of the Markov evolution equation in Lindblad form: 
j t Pt(B)=C(P t (B)), C{B) = -G*B-BG+ l -Y^L* l BL l: G = %H +\Y^L* l L l . 

This research was partially supported by the Russian Foundation for Basic Research under grant No. 95- 
01-00784 
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